Power-law scaling in dimension-to-biomass 
relationship of fish schools 
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Abstract 

Motivated by the finding that there is some biological universality in the relationship 
between school geometry and school biomass of various pelagic fishes in various 
conditions, I here establish a scaling law for school dimensions: the school diameter 
increases as a power-law function of school biomass. The power-law exponent is 
extracted through the data collapse, and is close to 3/5. This value of the exponent 
implies that the mean packing density decreases as the school biomass increases, 
and the packing structure displays a mass- fractal dimension of 5/3. By exploiting 
an analogy between school geometry and polymer chain statistics, I examine the 
behavioral algorithm governing the swollen conformation of large-sized schools of 
pelagics, and I explain the value of the exponent. 

Key words: power-law scaling, data collapse, pelagic fish, school size, geometry 



1 Introduction 



1.1 Brief history of studies on the packing geometry offish schools 



Animal packing in social aggregations is of fundame ntal interest in ecolog 
their conformations have been extensively studied (jParil Il927l: ISvmons 



pkubo and Chiang . 119741 
1980: .Partridge. . 198 0". 'l982 
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/, and 



1^ 



Graves, 1976; Pit cher and PartridgeL Il97f 
.Dill et alL 11981c lAoki and Inagakil . ll988l : |Misun( 



1993al : Parrish and Hamner (Eds.)l . 1997 ) . Pioneer tank observations of pelagic 



fishes (herring, sprat and mackerel) were conducted bv iParn (|l927l), and the 
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"persistently uniform density" of a school was noted. From experiments on 
schools of saithe Pollachius virens, herring Clupea ha renqus and cod Gadus 
morhu a cruising in a ten meters circular gantry tank, Pitcher and Partrideel 
validated that all the fish in a school occupy a volume proportional 
to Nb^, where is the number of fish and b is the mean nearest- neighbor 
distanc e (approximate to the mean fish-body length) in the school. Misundl 
(19933) reported from field observations that the number densities of herring 
schools are almost independent of the dimensional size of school but are an 
order of magnitude lower than the packing densities observed when schooling 
in small tanks. 



Taut i and c olleagues ( Tauti and Mivosi , 19291 : Tauti and Hudinol . lT92i l iTauti and YasudaL 
1929Lll93dir933a> .bl had experimentally shown that a fish school can be viewed 
as a group of inorganic particles and treated with the methods of physics. 
Brederl ()l954 ) explicitly applied physical equations to such fish schools. Re- 
cently, a number of theoretical and numerical models of schooling and flock- 
i ng behavior have been s t udied by biologi s ts, raathematicians and p hysicists 



( Okubo and Levinl . 120011 : ICamazine et al.L l200ll : IVicsek fEd.l 120011 ). As for 



the animal group geometry, many models for social aggregations, however, pre- 
dict that densities increase as th e group size (in number) increases [overviews 

of suc h models are given in Re fs. ( Warburton and Lazarusl . 1991 ; Beecham and Farnsworth , 
19991 )]. iMorilner et all (|2003h mathematically tackled this problem of constant 
density and revealed the condition for a well-spaced group, i.e. what class of 
mutual interactions results in a relatively constant individual distance in the 
interior of the aggregate. 



More recently, by means of underwater acoustics, the school sizes (in num- 
ber or biomass of fish) have been quantitatively measured for different values 
of the dimensional size of schools in the wild. Precise data on conformations 
of large-sized schools of pelegic fishes became available rapidly, which were 
extremely helpful in elucidating a certain geometric law, i.e. the relation be- 
tween dimensional and biomass sizes of pelagic fish schools, bringii ig about 
some i mportant cha nges i n our vi ewpoint s. Misund and colleagu es (iMisun^ , 
1990: Misund et al.l.J992:.Misund.. 1993b: iMisund et al.l . Il99,4 Il996. : .Coetzee . 
2OOOI: iMisund and Coetzed . I2OOOI : iMisund et al.L l2003lH ound that the power- 
law scaling in dimension-to-biomass relationship exists, and is robust across 
a broad range of pelagic species as well as across diverse environments. They 
demonstrated that if the biomass iV in a school is, say, doubled, the cross- 
sectional area of a school is increased by a factor 2^^^ 



I.e. 



cross-sectional area oc (biomass) 



2v 



and that the exponent u looks universal, reading 0.5 (from the field data u 
ranges from 0.415 to 0.77). The geometric law they found implies that the 
mean density of a school scales as N^~^'^ in three dimensions of space and the 
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conformation of social aggregations swells (i.e. v > 1/3). Such a relationship as 
Eq.(l) with u = 0.5 has been utilized for the dimensions to biomass conversion 
(e.g. transforming the sch ool diameter R to the number iV of fish) in analyz ing 
school-size distributions (|Andersonl . Il98l[ iNiwaL ll996bL ll998L 120031 l2nn4ah . 



It might come as a surprise that packing densities decrease as the group size 
increases, contrary to previous observations and predictions. Laboratory obser- 
vations for school geometry and internal structure have been made exclusively 
in small tanks, which generally show the constant density of fish that the school 
volume is proportional to the number of individuals. In all set-ups, the factors 
resulting in homogeneous, cohesive school may be especi ally pronounced, and 
i t ther efore seems that the quantified structure is skewed. iPitcher and Parrish 
f)l993t) claimed that homogeneity in schools has been over-emphasized. In-situ 
observations of herring and sprat schools with a high resolution sona r revealed 
that the packing structure within the schools is rather heterogeneous (jCushine , 



19771) . This has been confirmed by measurements of free-swimming schools us- 



ing photography and high-resolution echo integration, which showed that the 
packing density dist ribution in capeliri Mallotus villosu s and clupeoid schools 
varies considerably ( Freon et al. , 19921 : Misundl . 1993a ) . Regions of high den- 
sity are u sually found within the schools, and even empty vacuoles have been 
recorded. Misund and Floen ([l993i) observed by repeated echo integration that 
there were large variations in internal packing density of herring schools (i.e. 
high-density regions or empty lacunas within a school), and that the packing 
density structure was quasi-stationary. Besides artificial environments in small 
tanks, the discrepancy in former observations could ha ve been caused by too 
small numbers of fish in the schools. For instance, in Pitcher and Partridge! 
f)l979>) . N takes a few tens of fish. Since such a geometric law above is always 
defined only in a certain limit ( de Gennei . 19791 ). the scaling in the dimension- 
to-biomass relationship is expected to hold for large-sized schools of pelagics. 



In this paper, the exponent z/ is estimated according to the est ablished uni- 
versal scaling law in the school-size distribution of pelagic fishes fjNiwal . I2003L 



2004a|): choose the suitable value of v to achieve the best data-collapse on the 



size distributions in terms of the school dimension. Notice that the dimension- 
to-biomass relationship is a property of the single school and the scaling expo- 
nent V is determined by the behavioral algorithm of fish schooling at individual 
level, while the scaling in the school-size distribution emerges from the inter- 
school interactions at population level (i.e. a global property of the interacting 
school system). 
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Fig. 1. Vertical d imens i on (in meters) of Japanese sardine schools. Data from the 
acoustic survey of f sout heastern Hokkaido in the summer, 1982. The 

same data are available in Hara (il984[ ). Small schools of vertical dimensions in the 
first bin (< two meters) may be stray ones. 

1.2 Scaling in school-size distributions 



Animal group size is a focal issue in ecology that, in contrast to scaling, has 
introduced a single prefe rred size (i.e. optimal or compromise size) for any or- 
ganis m living in groups (^ PuUiani and Caracol 1984 : Higashi and Yamamura . 
19931 : iNiwa. 1996b; Hoarc e t all 120041) .1^ gure 1 shows an example of the his- 
togram of school dimensions of Japanese sardine Sardinops melanostictus. A 



peak frequency and a right skew are typical of pelagic fishes (|Andersonl . Il981 



Niwal . Il996bl . Il998l ). Their linear dimension, e.g. the vertical thickness, of the 
school falls into a certain range below a few tens of meters. One fish may 
not be the right atomic unit in schooling, since field observations suggest that 
no school exists under a certain minimal size. T his may cause b inned data 
of scho ol sizes to exhibit the fake peak frequency. lAnderson (1981) and Niwa 
f)l996b[ ) lave ignored slowly decaying (or fat-tailed) distributions, including 
scaling laws, in such data of school sizes short ranged with a fake peak fre- 
quency (the possibility of power laws was already presented in their models 
but not exploi ted). A tradi tional, widely used Gauss statistics says that, for 
the data from Hara ( 199Cll l shown in Fig.l, finding sardine schools ranging 
from 18 to 20 meters in vertical thickness should only oc cur about once every 
10^ detections of schools [for detail, consult iNiwal (j2004al )]. In other words, it is 
not the real world! Aquatic observations actually say that finding such schools 
occurs about once every 500 detections. The probability that such schools are 
found is 10^ times as large! 



The p ossibility of scaling in such distributions was found bv lBonabeau and Dagorn 
( 19951) . i^'at-tailed distributions have been found to quite generally character- 
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ize the size heterogeneity of pelagic fish scho ols ( Niwa . 19981 : Bonabeau et al. 



1998Lll999l ). Only lately. iNiwal ^200A 12004^ showed that the school-size dis 



tribution of pelagic fishes has the universal functional form of a power-law 
decay and a crossover towards the exponential distribution. 



1.3 Scaling and data collapse in ecology 



Scalin g laws have been found to characterize many patterns in e cological sys- 
tems |Azovskvl . l2nnn[lChave and Levinl . liooi iHallev et all 120041 ). There is an 
empirical rule about the relationship between the areas (L) of island and the 
numbers (S) of species on them: whether counting birds, butterfiies, plants or 
copepods in ponds, a tenfold smaller area contains roughly half the species. 
This can be fitted by a power-law function 



S{L) oc L' 



(2) 



with the exponen t z ~ 0.25 (|Prestonl . Il962l : iMacArthur and Wilsonl . Il967l : 
Harte et al.l . ll999l: iMavl . E999I ). The species-area relationship is widely used in 
rough assessments of likely future rates of species extinction (because observed 
rates of tropical deforestation are equated to loss of habitat area, we can as- 
sess the annual production of species committed to extinction). So, a mo re se- 
cure u nderstanding of s uch re lationships has great practic al importance (iMavl . 



1999 h iBanavar et al.1 (Il99 



I972I : iBinder and Heermam 



f appl ied finite-size scaling (jPisher and Barbei 
1988| ) to test the hypothesis of scalin g invariance 



i n the species-area relationship resting on a model proposed by iHarte et al. 



(1993) :br the species abundance distribution. Following the conventional no- 
tation in the physics literature, the system sizes are hereafter denoted by L. 



The finite-size scaling (FSS) hypothesis assumes that the fraction of objects 
with size n for a finite system of size L is written, with a scaling exponent A, 
as 

P{n; L) = L-^F (n/L^) , (3) 

where F is a universal scaling function. P{n; L) reads, for instance, the species- 
abundance distribution, which defines the probability that any given species on 
a census patch of area L has n individuals, where P{n; L)S{L) is the expected 
number of species with n individuals. is a measure of the width of the 
probability distribution P{n;L), e.g. the mean or the standard deviation of 
the distribution. Equation (3) expresses the principle that the behavior of 
the finite system, i.e. P(n;L), is controlled by the ratio n/L^. We test the 
FSS hypothesis by data collapsing: when y = L'^P^n; L) is plotted versus 
X = n/L^, if the distributions across different L's fall on a single curve (e.g. 
collapsing distinct abundance distributions measured in different areas and 
with different total numbers of individuals onto one scaling curve), then we 
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should be able to identify a scaling function F{x). In order to determine the 
exponent A, we try a best-fitting procedure such that the family of curves 
P(n; L) collapses onto a single curve as well as possible. We then obtain other 
scaling exponents (e.g. z of the specie s-area relationship) resting on the scaling 
relation between scaling exponents (|Goldenfeldl . Il992h . By noting that the 
total number of individuals of all species is equal to J^nT^Pi^T'i L)S{L), the 
species-area relationship leads to a scaling relation 



A + z = l, 



(4) 



assumi ng that the total nurn ber of individuals of all species is proportional to 
area L ( Banavar et al. . 19991 ). 



The universal scaling is fundamental to data collapse. The power-law scaling 
emerges as ubiquitous properties in ecology. In statistical analysis in ecology, 
data collapsing across different environments or species has been, however, ob- 
served only lately in dynamics of breeding bir d populations ([Ke itt an d Stanlev , 
19981 : Keitt et al. . 2002) . food web struct ure (ICamacho et al. . 12002). microbial 
body-mass spectra in marine .nn. v.t.^jRi..lHo.t„ll knnA school- size dis- 
tributions of pelagic fishes (|Niwal . l200il2004al ). and in the context of ecological 
economics, capture fisheries productions of countries ( Niwal 2004bl ). 



2 Scaling in dimension-to-biomass relationship 



2. 1 Estimation of the exponent v through data collapse 



In order to understand the geometric properties of school configurations, the 
school biomass must be measured for different values of dimensional size R, 
and we must compare them. Misund and colleagues extensively performed the 
simultaneous observations of the values of R and for pelagics by the follow- 
ing two methods: (i) two-dimensional (cross-sectional) acoustic measurements 
(unit in square rneters) and su bsequent purse seine capture (unit in tonnes) of 
schools (|Misundl . I 199o[ ligbsbl ). and (ii) acoustic measur ements of the three- 
dimensiona l structures and backscattered echo e nergy (iMisund et a l I Il99i 
199,4 Il99fil : ICoetzeel . I2OOOI: iMisund and Coetzed . l2000[ iMisund et allbOO.^ . 
The surveys were conducted on stocks of anchovy Engraulis capensis, her- 
ring Clupea harengus, horse mackerel Trachurus trachurus capensis, mack- 
erel Scomber scombrus, pilchard Sardinops ocellatus, round herring Etrumeus 
whiteheadi, saithe Pollachius virens, sardine Sardinops sagax, and sprat Sprat- 
tus sprattus, in different seasons and geographic regions (the Barents Sea, the 
North Sea, the Norwegian Sea, the northeastern Atlantic and off Namibia, 
and off the coast of South Africa). They found that there is some biological 
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universality in the dimensions-to-biomass relationships, Eq.(l). Their finding 
implies that the radius R of the school scales as 



R = (constant) x A^'^ 



(5) 



in a statistical sense. N, denoting the school biomass, is hereafter defined by 
the number of fish in a school. 



Now 1 propose another way to establish the scaling law for school dimensions, 
Eg. (5) , by applying universal scaling law in fish school-bio mass distributions 
( Niwal . l2nn.il2nn4at ). The biomass distributions W{N) follow a power law with 
exponent /3 = 1 up to a cut-off size {N)p, 



W{N) = N-^P {N/{N)i 



(6) 



where P{x) is a crossover scaling function with a strong drop for x > 1, and the 
cut-off size (crossover size from power-law to exponential decay) is calculated 
from the biomass histogram data {{Ni, Wi)\i = 1,2, .. .}, 



where histogram bins are chosen with width AA^. 



(7) 



The dimension data of fish schools are binned with width AR, giving the set of 
frequencies | (^Ri, W^''^^^ i = 1, 2, . . .|. From Eqs.(5) and (6), the distribution 
of geometric dimensions of fish schools is represented as 

W^^'^R) = R-^P^""^ {R/{R)p) , (8) 



where 



and 



{R)p = (AT) 



p(G), 



X 



Pix^h 



(9) 



(10) 



Therefore, the school-dimension distribution follows a power law with the same 
exponent "—1" as the school-biomass distribution. The power-law distribution 
W'-'^\R) is truncated at a cut-off size {R)p, which is also calculated from 
histogram data of school geometric dimensions. 



{R)i 



EiRf^'^W.^^'^AR 



The following normalizations are adopted for the scaling function P'^^\x) and 
the histogram data of geometric dimensions of fish schools, because the cut- 
off size (N) p is proportional to the total number of fish in the school system 
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Niwal (|2nn.il2nn4ah : see also Eq.(23)]: 

Jo 



and 



respectively. 



j:Ry''wrAR={Ryj 



l/u 



:i2) 



(13) 



Fat-tailed school-size distributions are necessarily truncated because the pop- 
ulation is finite. Since the size {R)p depends on the exponent u following 
Eq.(ll), so that we can determine the value of u by evaluating the location of 
the cut-off in the power-law distribution W^'~^\R). The ordinary least squares 
regressions raight, however, lead to a "w rong" estimation of the exponent 
( Niwal . ll998l:lBonabeau et al.l . ll998L IToOfll ). '. make use of the data collapse to 
extract the "right" exponent. From Eqs.(8) and (12), when y = W^'^^{R)p is 
plotted against x = R/ (R) p with correct parameter z/, all the empirical data 
should collapse onto each other. The power-law exponent of school-dimension 
distributions, u, is then evaluated through data collapse. Let us search for the 
value of z/ that places all the points most accurately on a single curve. We use 
a set of histogram data of vertical dimension of J apanese sar dine Sardinops 
melanostictus schools, from 22 acoustic s urveys bv iHara (|l99nl ) off southeast- 
ern Hokkaido for seven years, 1981-1987. Hara ( 1986t ) reported that Japanese 
sardine migrate as a huge-sized school in number from a few hundreds of thou- 
sands to a few million of fish. To obtain the best data collapse, the x-axis is 
divided into bins (Fig. 2a), and for each bin two-dimensional variance 



(14) 



is calculated, where and ay denote the standard deviation of the mean fi^. 
and fly, respectively. The parameter u is then estimated at value that minimize 
the mean of two-dimensional variance for the bins (Fig. 2b). The mean of two- 
dimensional variance, e , is a measure to determine the goodness of collapse 
( Lillo et al. . 20021 20031 : Niwa . 2004a ). A good data collapse can be obtained 
by using the value u ~ 0.6. The resulting plot of empirical school data is 
shown in Fig. 2a. Experimentally fitting the parameter u to achieve a good data 
collapse, "3/5", is reminiscent of the Flory value of the exponent in a power- 
law dependence of the coil radius of a pol ymer chain (in three-di r nensio nal 
solutions) on the degree of polymerization ( Florv . 19531 de Gennei . 19791 ). 



Here w e see that t he po wer-la w regi n ae of th e distribution is too short, which 
misled Anderson (1981) and Niwa ( 1996bl ) into overlooking the power-law 
distributions of school sizes. Notice that many power laws that appeared 
i n the ecolog y literature span less than two orders of magnitude of scale 
( Hallev et al. . 20041 . The power law range of too few scales is not unique 
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Fig. 2. Data collapse to extract the exponent v. (a) Scaled distribution of geometric 
dimensions of sardine schools, y = W'-'^\R)p is plotted versus x = R/{R)p with 
I' = 0.6 on double-logarithmic scale. The bins are chosen equally spaced on a log- 
arithmic scale as X G |^lo^i+'^/^, 10^^+^^+^)/^) with k = 1,2, ...,6. The rectangle 
in gray reads the interval Hy ± ay, i.e. the error ay on the mean fiy (indica ted b' 



the s lit) for each bin. The solid line is a prediction of the mean- field theory ( Niw; 



I 



2003). (b) The mean of two-dimensional variance, e, versus the power-law exponent 



v. Al though e sh ows noisy fluctuations, it takes a minimum around v = 0.6. Data 
from iHaral (jl99flh are analyzed. 

to ecology; the larges t numbers of power laws reported in the physical science 



are for small ranges (jHamburger et al.L |1 



2.2 Retest of the FSS in school-biomass distributions 



The acoustic-survey data are converted into a school-biomass histogram as 
follows 

W{N)AN oc W^''\R)R^/''-^AR. (15) 

We now crosscheck the value of u through finite-size scaling (FSS) analysis 
of school-biomass distribution ( Niwal . 2004a ). Since the finite population size 



causes the truncation of power-law distribution W{N) oc A^~'^, there is a well- 
defined quantity 

L= \ ' , (16) 
EiNj'WiAN 

which depends on the system population size. In order to characterize the 
finite size effects, FSS hypothesis is used: the distribution function depends 
on N only through the ratio N/L^, 

W{N; L)dN = L-^F [n/L^) d (n/L^) , (17) 
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where F{x) is a universal function independent of fish population size. The 
prefactor is required to ensure the normalization 

J2K^W,AN = 1, (18) 

i 

where P{N) [= N^W{N)] now represents the probability distribution of 
school-biomass sizes. From the FSS hypothesis, it is expected that when 
]yijA+B jg piQ|;|;g(j against N/L^ with correct parameters A and B all the 
data collapse onto a single curve. The power-law exponent of biomass distri- 
butions, P, is then evaluated through FSS analysis. The value of B/A is the 
estimate of the power-law exponent 

(19) 

Let us analyze a set of 22 acoustic-survey data of sardine schools (same as 
Fig. 2) converted into biomass histograms by using Eq.(15) with u = 3/5. In 
a simultaneous best-fitting procedure (Fig. 3), a good data collapse can be 
obtained by using the values A ^ 1 and B ^ 1. The power-law exponent 
derived from the FSS collapse is /? ~ 1. The resulting plot is shown in Fig. 3a. 
The school-biomass distribution follows a power-law decay with exponent — 1, 
and is truncated at the cut-off size of Eq.(7). The FSS collapse confirms the 
scaling laws for school sizes, Eq.(5) with u = 3/5 and Eq.(6) with P = 1. 

Finite-size scaling techniques have been applied to ecology and shown to 
predict scaling r elations betwee r i scal i ng exponents in relative a bundance of 



species [Eq.(4)] fiBanavar et al. . 19991: Mi and Gold enfeldL 2001 ). in dynam- 
ics of breeding bird populatio i is (iKeitt et al.L 12002), in biomass-size di stribu- 
tions o f seston (jRinaldo et al. . 2002[ ). and in exploitation of fish stocks (Niwa, 
2004bh . In school-size statistics, we expect to have a scaling relation. We now 



choose the normalization 



Y^N,W,AN = <i>, (20) 

i 

where $ denotes the total fish population in the school system (J2i WiAN 
gives the total number of schools). Since Eq.(16) implies oc with a 
scaling exponent 7, the FSS relation for the school-biomass distribution is 
written as 

W{N) = N-f^P^N/^^). (21) 
Accordingly, the normalization of Eq.(20) yields the scaling relation 

2 -(3 = 1/7. (22) 

The best-fitting value in data collapse, /3 ~ 1, gives 7 1. As a consequence, 
the location of the cut-off in the power-law distribution of school sizes simply 
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1.03 




Fig. 3. FSS analysis of school-size data, (a) FSS plot of the biomass distribution on 
double-logarithmic scale. Dimension data of sardine schools (same as Fig. 2) are con- 
verted by Eq.(15) with v = 3/5. Here y = WL^~^^ is plotted versus x = N/L"^ with 
A = B = 1. Two-dimensional variances [same as Eq.(14)] are calculated for bins 
chosen equally spaced on a logarithmic scale as x G ^W'^+^'^'^-^y^ ^ ]^Q-i+(fc+o.5)/3^ 
with k = 0,1, ... ,5. The rectangle in gray is same as Fig. 2. The solid line is a pre- 
diction of the mean- field theory (Niwa. 21)0.^) . (b) The region of the AB-p\ane in 
which the minimum of the mean of two-dimensional variance exists. A measure of 
data collapse for scaling, i.e. the mean of two-dimensional variance, e, takes a min- 
imum Emin for the right choice of {A,B). The minimum is found with the precision 
i.e. width of the minimum, Ae = 10~^ in black region (Ae/emin ^ 3.12 x 10"'^). The 
values of the parameters lie in the intervals A = 0.999±0.008 and B = 1.003±0.024, 
and therefore /3 = 1.004 ± 0.032. 

reads 

{N)p oc $, (23) 
which is verified by means of extensive numerical simulations ( Niwal I2C 



3 Behavioral algorithm of fish schooling 



3.1 Gaussian model 



Let us now investigate cohesive motion of schools in a three-dimensional space 
from the viewpoint of the behavioral algorithms which govern their formation 
and dynamics: attraction between neighbors maintains cohesion of the school; 
a tendency to align with neighbors produces collective motion of the school. 
The minimal model of cohesion is a linkage of neighbors consisting of har- 
monic spring, because the inter-fish distance follows a Gaussian distribution 
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( Dill et al. , 1981 ; Partride^ . 19821 : Niwa , 1996a ) . Using the relative coordinates 
to the center of school, the equation of motion of fish (as noisy self-propelled 
particles) in a large school of size N is wr i tten by a on e-body approximation 
as the following Langevin equation (Niwa, 1994 . 1996a| l: 



(24) 



providing that individuals are sufficiently sensitive to behavior of their neigh- 
bors, where f{r) is the cohesive force; J denotes the strength of alignment; and 
77 is a random perturbation of the velocity with strength and a 5-correlated 
time dependence. 



3.2 Self- consistent calculation of of cohesive force 



The attractive interaction force acting on one body due to a system of bodies 
is the neighbor- joining harmonic spring. One essential approximation is to 
replace the many-fish problem by the problem of solving the motion of one 
fish in a certain self-consistent field. The cohesive force is then written in the 
following form by the one-body approximation: 



fir) 



(25) 



where kN'"" is the effective spring constant for the springs strung out from 
the center of the school to a fish: the number of consecutive springs joining the 
fish to the center of the school via other companions is proportional to A^". 
The total "elastic energy" of inter-fish bonds in the school depends linearly 
on number of individuals in the school. The overall elastic energy after 
integration over a sphere of radius R (denoting the average radius of the 
school) results in E^i oc R^N'"", which is derived from the ideal harmonic 
spring, Eq.(25). Therefore, the exponent a will be determined self-consistently. 



(26) 



The solution of th e Langevin equation (24) with Eq.(25) takes the following 
asymptotic forms (|Horil . Il977l ): 



o"», 



{dr/dty 



(27) 



2 

k ' 



f28) 



where the root-mean-square y (r^) gives the average radius of the school, R. 
From Eq.(26) together with Eq.(28), the self-consistent value of the exponent 
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a is obtained: 

a = 1. 

As a consequence, the self-consistent cohesive force reads 



fir) 



(29) 



(30) 



where h denotes the effective inter-fish distance (a constant independent of 
A^): 6^ = al/k. The average radius of the school is then given by hN^/^. Note 
that Eq.(27) is an ex ample of a m ore general principle called the fluctuation- 
dissipation theorem (lKubol . ll96fit) . 



Notice that the above model has been developed in the absence of the ex- 
cluded volume interactions. Social cohesion by harmonic spring interactions 
between neighboring fish reminds us of the classical picture of a polymer 
chain based on the bead-spring model (|Rousel . lT95l . The exponent v close to 
the value 3/5 may then be understood by taking the excluded volume effect 
f Poi and Edwarda . 1986t ). Such a school without excluded volume effect may 
be called the "Gaussian" school in line with polymer physics. 



3.3 The excluded volume effect 



From Eq.(24), the probab ility of the school radius being between r and r + dr 
is given by the following (|Hori Il977l l : 



^o(r) =47rr2 f-^^jexp (--^ 1 (31) 



in a stationary state in three dimensions (i.e. the position vector r follows a 
Gaussian distribution). The Gaussian school model considered above permits 
fish to occupy the same region in space. Of course this is a physical impossibil- 
ity since each fish possesses its own finite volume. Therefore, in the "excluded 
volume" school, there are a number of Gaussian school configurations which 
are disallowed due to the steric effect. Let p{r) be the probability that a Gaus- 
sian school configuration, as counted in Eq.(31), is also allowable under the 
excluded volume condition. We now calculate the probability that no overlaps 
occur when we place N fish within a region of volu me (^ r^), which will lead 
to an estimation for p(r). The approach is due to IPoil ( 19961 ). Letting w be 



the volume which is effectively excluded to one fish by the presence of another 
{w < 6^), the probability that one particular fish will not overlap with another 
is then given by {1 — w/r^). Since there are N{N — l)/2 possible combinations 
of pairs, the probability that no overlap occurs in all of these combinations is 
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given by 



p(r) = (1 - u;/r3)^(^-^)/2 ^ (33) 



where 3> w and ^ 1. Therefore, the probabihty distribution of the 
school radius r can then be written as 

^(r) = ^'oWpW ocr^exp I - ^-^j (33) 
for the excluded volume school. 

Both "^oir) and \E'(r) have a maximum at certain values of r. Let us estimate 
the average size of the school radius in each model by calculating the positions 
of these maxima. The maximum of "^oi^r) occurs at Rq = (26^A^/3)^/^. The 
maximum of \l'(r) occurs at R, which satisfies the following equation obtained 
by differentiating the logarithm of Eq.(33): 



R\^ fR\^ 9VQw^/N ^^^^ 
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If iV 3> 1, the second term on the left-hand side of Eq.(34) can be neglected 
and hence 

1/5 



i?^i?o(^) ocN^/'. (35) 

The exponent 3/5 is the very value extracted through the data collapse. 

The scaling in the relationship between geometric dimensions and biomass 
of pelagic f i sh sc hools is analogous to that developed in polymer physics 
fide Genne.4 Il979h . The result suggests that the dynamics are unexpectedly 
"simple" and depend primarily on common cohesive motion in these animate 
and inanimate systems. Though the above is a very rough theory of school 
conformation with the excluded volume effect, the overall statistical proper- 
ties do not depend on the details of the model, which is a consequence of 
universality ( Stanlev . 1995[ ). 



3.4 Exploiting the analogy to a polymer chain model 



As indicated in the Int roduction, large int e rnal v ariations in packing density 
occur within a school. Freon and MisuridI ( 1999| ) pointed out that a source 
of substantial variation in internal school structure is the formation of sub - 
groups. Such subgroups have been observed in sait he schoo l s (IPartridgel 1981 ). 
and in schools of minnow P hoxinus phoxinus (Pitcher, 1973| ) and herring 
( Pitcher and PartridgeL 19791 ). Relatively independent movements of such clus- 
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ters of individuals can open up empty spaces and cause large variation in school 
volume. 



The pol ymer-chain analo gue of the subgroup in school conformation is the 
"blob" (jde Genne.4 1197^): the polymer chain behaves as a series of blobs. 
Based on the blob concept, for small value of the excluded volume w, a sub- 
group within a school, with a number of Ab of fish, must be nearly Gaussian. 
We see this from Eq.(33) when we find no effect of w if wX]/'^ /b^ < 1. There 
is a certain value of Ag [~ [b^/w)"^], beyond which excluded volume effects be- 
come important. A single school will appear Gaussian at scales r < Tb where 



Tb ~ bX]/'^ ~ b'^/w, 



(36) 



while at scales r > Tb it will show excluded volume effects. According to 
the "blob" approach, an A^-sized school can be described as a cluster of N/ Ab 
subgroups. Inside the subgroup the core repulsion (by excluded volume effects) 
is a weak perturbation leading to a Gaussian state. The school radius is then 
written as i? ~ rB(A~"^A^)'^ in the native state. 



Notice the plasticity of school geometry and internal structure. The compact, 
dense packing, as observed in artificial environments like small tanks, may 
be caused by the very strong stress. When the school is confronted by dan- 
ger such as predators, the inter-fish distan ce decreases rapid ly and all the 
vacuoles within the school collapse quickly ( Freon et al. . 1992| ). This change 
may be closely related to the coil-glo bule transition in a polymer solution 
( Ptitsvn et Il968l : IdTGennej . Il97,'^ . 



In real schools the nature of the short range interaction is quite complicated 
like van der Waals-type inter-molecular forces. Accordingly, all fish within a 
school interact via the two-body potential deformed from a parabolic potential 
of harmonic spring [grouping forces have been expressed as th e grad ient of 
a potential function in previous model studies ( Morilner et al. . 2003| )]. The 
potential change, denoted by AUij, is given as the correction term to cohesive 
interaction potential between the i-th and the j-th fish: AUij will include 
steric effects and also may involve weak attractions. The correction factor to 
the probability distribution \I^(r) of the school radius r, in place of Eq.(32), is 
then given by the configuration integral 



exp I -J^AU,j/eu ) dxi ■ ■ ■ da; at, 

i<j 



(37) 



where Xi denotes the position vector of the i-th fish, and exp(— JA[/jj/e^) is 
integrated over the scho ol configuration space. We can estimate p(r) by means 
of the cluster expansion ()Maver and Maveii . Il94n[ lUhlenbeck and FordLE962l 
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Introducing a function of separation between 

Cij = exp (- JAf/,,/e,) - 1, (38) 

we integrate Eq.(37) in the Mayer cluster expansion over diagrams up to the 
second irreducible cluster, yielding 

^'W^^^p(-^-^^) (39) 
with the excluded volume parameter 

w = - I Ci2dri2, (40) 



and the other parameter defined by 

Wa = - yy C12C13C23 driadrig. (41) 

The actual configuration is decided by maximizing the distribution function 
^(r) = \l/o (r)p(r), and this maximization condition leads to the equation 
determining the school radius i?, 

\Ro) \Ro) 1666 [rqJ 16 63 • ^ ) 

Notice that, when the first irreducible cluster integral in Mayer's expansion 
is only taken into account, the configuration integral, Eq.(37), is evaluated 
by Eq.(32). In case that AUij consists of a hard core repulsion and a short- 
range weak attraction (i.e. at the distance of near-collision (r ~ w^^^) the 
attraction gradient between fish becomes steeper), Eq.(40) is estimated as 
w = wo{l — ee / €v) , expressing the dependen ce of w on the velo c ity-flu ctuation 
of fish, where Wq and ee are constants (jPoi and Edward.^ . Il986t ). During 
predation threats or laboratory observations, the value of excluded volume 
parameter w will change sign from positive to negative, provided that the 
velocity fluctuation is suppressed below a certain value ee- Decreasing e„ 
below ee, the dimensional size of the school becomes much smaller than that 
of a Gaussian school, as depicted in Fig.4. According to Eq.(42) the so called 
"expansion factor" R/Rq is determined not by w but by wN^^"^, and so if N is 
large only a small change in will cause a big change in dimensional size. For 
example, for school of 10^ fish, a small variation in will induce a dramatic 
change in the school radius. Equation (42) states that when —wN^^'^/b^ ^ 1, 
R/Rq is proportional to (— wA^^/^)"^/^ a^cl hence the solution is as follows: 

^-u;6W 
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Fig. 4. Swollen-dense packing transition in the school geometry [the solution 
of Eq.(42)]. Here R' = R/Rq is plotted as a function of a combined variable 
w' = (9\/6/16) wN^f^ /h^. Shown is the transition region \w\N^/'^ /h^ ^ 1. At = ee 
the excluded volume parameter w equals 0. The repulsive excluded volume effect 
then balances the attractive forces between fish, and the school behaves as a Gaus- 
sian school. The property of the swollen-dense packing transition is quite different 
in the two regions separated by the critical value w^/h^ = 0.0228. The transition 
between the the swollen conformation and the compact packing is gradual when 
w^/b^ > 0.0228 (dotted line), while it becomes extremely sharp and the transition 
is discontinuous when W/^/b^ < 0.0228 (gray line). The solid line shows the marginal 
transition behavior, and the crossover takes place when w' = —0.242 (indicated by 
the arrow). Notice that the solution of Eq.(42) is multivalued [depicte d by a sigmoid 
curve (white cutout)], showing a first-order discrete phase change (jPtitsvn et al.l . 
19681 : Ide 06x^^ . 1197,^ . when w^/b^ < 0.0228. 



which shows the closest packing of fish in a school, i.e. the constant density 
independent of the number of individuals. This change may be called the 
swollen-dense packing transition of schools. 



3.5 Configuration inside a swollen school 



The approach of treating fish schools as interacting particle systems naturally 
leads to the idea of applying successful methods of stat i stical physics to the 
description of moving together without a leader (jVicsekl . 120011 ). In statistical 
physics, the presence of non-trivial scaling is usually taken to mean that the 
dynamics are largely governed by simple geometric properties of the system 
and do not dep end strongly on detailed properties of the system components 

( WiisonL irr^ 



As for the school geometry, the basic units are the effective inter-fish length 
b and the number of fish, N. We assemble A neighboring fish in a school into 
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Nl A groups, e . g. in t he sense of the Voronoi-Dirichlet diagram (|Dirichletl . ll85n[ 



Voronoi . IQOSl 19091 ). The length between centers of neighboring subgroups is 



defined by X'^h with an exponent v. After such a scale change as 

A^ ^ A-^A^, h X%, (44) 

the macroscopic quantities which determine the overall properties of the school 
satisfy 

giX-'N,X%)=X^giN,b), (45) 

where x is an exponent which depends on the quantity under consideration. 
For this to hold true for arbitrary A, the function g{N, b) must take the scaling 
form 

^(A^, b) = N-''g^'''\bN^), (46) 

where a new function g^^'^^ is introduced. Th is is the Widom-Kadanoff scaling 
law in critical phenomena ( Goldenfeld . 1992f ). For example, the average radius 



of the school, R, should be unaltered under the transformation of Eq.(44): 

R{N,b) = R{X^^N,X%) = R^''\bN''). (47) 
Because of Eq.(5), R^^'^\x) = x and we must have 

u = u, (48) 

where i/ = 3/5 in three dimensions. 

As an example of the application of the above scaling law, let us consider the 
pair correlation function inside a school, where fish movement take place in a 
three-dimensional space. A pair correlation function p2{r) is defined as follows. 
We pick one fish at random in the school, and we place it at origin. Then we 
ask, what is the number density of other fish at a (directional) distance r 
from the first, and we average the result over all choices of the first fish (and 
of directions). The function p2{r) has an integral which is just the total number 

of fish per school, J p2{r)dr = N. Hence, from dimensional analysis, we can 

write 

P2{r) = R-^g{N,r/b). (49) 

Under the transformation that neighboring A fish are grouped to form one 
subunit in a school, P2{r) will be reduced by 1/A, since the pair correlation 
function is proportional to the number density of fish. Therefore 

g{X-^N, r/X^b) = X'^g{N, r/b). (50) 

Then, the function p2{r) obeys a simple scaling rule 

P2{r) = NR-^pt\r/R), (51) 
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where p2 is a dimensionless universal function, and A^i?~^ is the dimensional 
factor. Focusing on the limit r <^ R, we can reach the form of p2{f') by a 
simple argument. In a sphere of radius r we have a certain number of fish, N, 
related to r by the excluded volume exponent: N^^^b ~ r. The function P2{r) 
scales like the density of fish in the sphere, 

N 1 

P2W^^^;:i75^ (for r<i?), (52) 
which gives P2'^\x) ^ x^^^^. 

The possible observable property, Eq.(51), tells us that if we were to measure 
P2{t) for fish schools of different size N, there would be superposition of the 
curves obtained by plotting p2{r)E? /N against r / R. Thus this kind of scaling 
relation may be verified experimentally on large schools of pelagic fish. 



4 Discussion 



The power-law scaling generally exists in dimension-to-biomass relationship of 
pelagic fish schools in nature. Here I have estimated the power-law exponent v 
of the geometric relation, based on the dependence of the power-law regime of 
the school-dimension distribution W^'~^\R) on z/, i.e. Eq.(ll). We have tested 
whether the distributions W^'^\R) are self-similar (i.e. exhibit scaling), rely- 
ing on both the power law for the dimension-to-biomass relationship of Eq. (5) 
and the FSS relation in the power- law school- biomass distribution of Eq.(6). 
Plotting the scaled histogram-data from the 22 in-situ observations, we have 
found that the 22 curves do indeed collapse onto each other (Fig.2a), suggest- 
ing that Vr^°'(i?) follows a universal functional form [Eq.(8)]. We have ex- 
tracted the power-law exponent of dimension-to-biomass scaling relationship 
via a minimization of a measure to quantify the nature of scaling collapse. 
Next I have explained the value of the exponent proposing the Gaussian 
school model for the fish with excluded volume w. By exploiting the analogy 
between fish-school and polymer conformations, we have examined the behav- 
ioral algorithm governing the swollen conformation of large-sized schools. We 
have seen that the exponent is modified strongly by the steric effect. What is 
universal in Eq.(5) is the exponent v ~ 0.6: it is independent of species as well 
as environmental conditions, and the same for all schools. The constant that 
multiplies N'^ in Eq.(5) is non-universal [~ (w6^)^/^ with effective inter-fish 
distance 6, predicted from Eq.(35)], and depends on the details of interac- 
tions between f ish which may vary with species and environmental conditions 
( Morganl . lT988h . 



To understand the scaling law for school geometry, Eq.(5) with z/ = 3/5, it is 
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essential to see what the value of u represents. Obviously, for a regular object 
embedded into a d- dimensional Euclidean space, Eq.(5) would have the form 
N{R) ~ expressing the fact that the volume of a (i-dimensional object 
grows with its linear size R as W^. Contrary to an integer dimensionality, it 
has been shown that the packing structure within the schools is characterized 
by a non-integer (i.e. fractal) dimensionality of 1/z/ ^ 1.7: the number of fish 
in a school of radius R scales as N(R) ~ R^^'^. During the last decades of 
the twentieth century it has widely been recognized by researchers working in 
diverse areas of science that many of the structures commonly observed possess 
a rather special kind of geometrical complexity; t he name "fractal" w as coined 
by Benoit Mandelbrot for these complex shape s (lMandelbra3. 19821) . Objects 



of biological origin are many times fractal-like (|Vicsek f Ed.! 120011 ). 



FSS plots of data lead to possibly important and practical insights in ecology. 
In the context of the fisheries mission, the demonstrated geometric relation 
between dimensions and biomass of pelagic fish schools, when applied to map- 
ping pelagic schooling fish, will largely improve the precision in the fish stock 
assessment. Because of the linear dependence of the fish population on {N)p, 
as given by Eq.(23), the fish stock abundance can be inferred from an in- 
dex (i.e. cut-off size) that can be determined directly from observations [see 
Eq.(7)]. 

I am very grateful to Leah Edelstein-Keshet for helpful suggestions. 
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